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Abstract 

In this paper we present a new approach how to extract polarization-dependent exchange- 
correlation energy densities for two-dimensional systems from reference densities and energies of 
quantum dots provided by exact diagonalization. Compared with results from literature we find 
systematic corrections for all polarizations in the regime of high densities. 
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I. INTRODUCTION 



Density functional theory (DFT)^^ and spin-density functional theory (SDFT)^ are 
powerful techniques to investigate interacting electron systems. However, the results of 
these methods sensitively depend on the quality of the approximation of the exchange- 
correlation (XC) functional by the local (spin) density approximation (L(S)DA) and/or 
gradient corrections. In contrast to three-dimensional (3D) systems where a large number of 
parameterizations for XC energy densities and gradient corrections is available the situation 
is different for two-dimension al (2D) systems. The majority of all calculations for 2D (e.g. for 
QDs see Refs. ^j^liyyilJllljIlH) relies on the parameterization of the XC energy density 
from Tanatar and Ceperleyi^ (TC). They numerically calculated XC energy densities as a 
function of the carrier density (density parameter r^) and spin polarization ^ by applying 
variational Monte Carlo and fixed-node Green's function Monte Carlo methods. The results 
for = (no polarization) and C, = 1 (full polarization) and r^-values 1, 5, 10, 15, 20 etc. 
were fitted and interpolated using a Fade approximation 
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with parameters aj for ^ = and 1 respectively (see Tab. IV in Ref. ^3). In the limit of — 
the asymptotic behavior is characterized by ao 1 + bvs + O (r^/^^j whereby logarithmic 
contributions (see ring-diagram summation in Ref. [2^) such as r^lnrg are not taken into 
account. The construction of the dependence on the polarization ^ is analogous to the 3D 
case^ where the interpolation of the exchange (X) energy in Hartree-Fock approximation 
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Attaccalite, Moroni, Gori-Giorgi, and Bacheleti^ (AMGB) pursue a similar approach. 
They calculate the ground state (GS) energy of a 2DEG as a function of and ^ with fixed- 
node diffusion Monte Carlo methods which also take into account back-flow correlations. 
Their parameterization of the XC energy density 



Bxc^'^'ir,, 0/Ry = ^x(r., 0/Ry + ^c(r., 0/Ry 
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comprises the well-known X energy from the ring-diagram method^^ 
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and the correlation energy 



1) 4'H^.,0/Ry + 2ao(r.) + 2a^{r,)e + 2«2(r.)e', (5) 



with 



et\rs.O = ex(r.,0 " (l + + j^^^) ^x(r.,0) (6) 
being the Taylor expansion of the X energy with respect to ^ of the order six and higher. 
This kind of representation of the XC energy allows for the identification of the term 2ao{rs) 
with the correlation energy ec{rs,0) and of the term 2ai(rs) with the spin stiffness. The 
a;i(rs )-parameterizations are generalizations of the Perdew-Wang forn>2^ to 2D 
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The parameters which were fitted to the numerical results are summarized in Tab. II of Ref. 

Q. 

Another recent work on 2D XC energy densitie82i is based on a purely analytical approach. 
The interaction-strength-interpolation (ISI) connects the two limits of strongly and weakly 
interacting 2DEGs. The result for the XC energy density reads 
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The functions and parameters entering the interpolation are 
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Two observations can be made with respect to these different forms of the XC energy 
densities: 

1. To visuahze the different parameterizations we show in Fig. ^ the relative deviations, 
i.e. £^xc/^xc ~ Figs. H (a) and (b) compare the AMGB- and TC-XC energy density 
on different scales for r^. We find a good agreement (f^xc'^^/^xc ~ ^ ^ 0.008) which 
is not surprising as both expressions are based on the same method. On the other 
hand it is nevertheless remarkable because AMGB explicitly calculated polarizations 
< < 1 whereas TC interpolated over this regime. Larger deviations can be seen 
for ^xc/^xc ~ 1 in (c) and for ^xc/^xc*^^ — 1 in (d). There, the relative differences 
go up to 0.04, especially for full polarization. This deviation should have its origin in 
the different applied methods (ISI analytical, TC and AMGB numerical). 

2. Besides looking at the difference between the various XC energy densities we can also 
examine their performance when they are used as input in DFT. In a recent paper— the 
agreement between SDFT/LSDA (using the TC parameterization) and ED densities 
for QDs was investigated and found to be fairly good. However this has not to be true 
for all systems. For example in Fig. |21we compare GS densities of a QD with = 8, 
uj = 3meV, S = Sz = 4:, L = 0{h = l in this paper). Especially for the density 
maximum and minimum the difference between ED and SDFT/LSDA is considerably 
large. Further, this effect does not depend on the applied parameterization. 

These two findings are our motivation to check the accuracy of XC energy densities in 
2D by following a new approach. Our starting point are GS densities and energies from ED 
from which the XC energy densities are extracted in two steps: First we design an iterative 
scheme to calculate the self-consistent XC potentials of the KS equations which reproduce 
the exact densities. The second step deals with the extraction of XC energy densities from 
these XC potentials assuming L(S)DA. 

Our paper is organized as follows: We briefly introduce the QD Hamiltonian in Sec. II 
and the SDFT/LSDA in Sec. III. Then we describe how to extract XC potentials from GS 
densities (Sec. IV) and XC energy densities from XC potentials (Sec. V). Sec. VI contains 
the results for systems with arbitrary polarization. We end with a short summary of the 
most important conclusions in Sec. VII. 
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II. QUANTUM DOT HAMILTONIAN AND GROUND STATE DENSITIES 



We consider a two-dimensional QD with an axially symmetric parabolic confinement 
potential of strength ujq. As we are especially interested in GS configurations with vanishing 
angular momentum L = a magnetic field can be omitted. The Hamiltonian for particles 
in real-space representation [with r = (x,?/), p = {pxjVy)] reads: 

Here m* is the effective mass, e is the electron charge, and e is the screening constant of the 
host semiconductor. 

In order to get the reference densities and energies for the investigation of XC energy 
densities we apply ED techniques which provide results of high accuracy.— Please note that 
the (spin-) density of all eigenstates of the angular momentum operator are functions of radius 
r but not of the angle if. Therefore, the relevant quantities from ED are the spin-densities 
n„{r) of the GS and its energy Eq. 



III. SPIN-DENSITY FUNCTIONAL THEORY AND LOCAL (SPIN-)DENSITY 
APPROXIMATION 

In this Section we sketch the basics of (S)DFT and L(S)DA. The DFT formalism was 
originally established by Hohenberg, Kohn, and Sham^"^ and generalized to spin-polarized 
systems^ by including the coupling of the polarization to an applied magnetic field. Ac- 
cordingly, the Hohenberg-Kohn (HK) theorem has to be modified with respect to the spin 
degrees of freedont^. For this case, it states that two different non-degenerate ground-state 
wavefunctions and always yield different combinations (no-(r)) ^ (n^(r)) of spin 
densities. This is sufficient to establish a functional of the total energy with the usual 
functional properties 

EvAn.] = FhkK] + E / drK(r)n.(r) (16) 
and the universal HK functional 

FhkK] = (^K]|r + w\nna]). (17) 

Thus, even in the limit of vanishing magnetic fields the SDFT scheme can yield a spin- 
polarized ground state for even electron numbers due to Hund's rule.^^^ 
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For practical purposes the variational scheme has to be mapped on the Kohn-Sham (KS) 
system, i.e. an effective single-particle system with the same GS densities as the interacting 
system. The spin-degree of freedom is considered in the KS equations^ by assuming the 
total spin Sz in z-direction to be a good quantum number 

r 



2m' 



+ K(r) + 4^ / dr'^^ + Vxc,.(K], r)| (p,,.(r) = (18) 



with the spin a = ± in 2;-direction and the KS energies 6i^„ < £2,0- < ■■• • For a system 
containing N particles the spin densities are given by 

^<x(r) = Il7i,a|<^j,<T(r)|^ (19) 
j 

with 7j_o- being occupation numbers of the KS levels in the ground state {J2j = and 
A'^.i. + = A^). Then the GS density and polarization are 

n(r) = ?2+(r) + ri_(r) (20) 
e r = ^ \ , ^ ■ 21 



The XC potentials 



W(K]>r) = ^-7;^ (22) 



are defined as functional derivatives of the XC energy functional 

ExcK] = FukK] - / dr I dr'^^^^^ - TsK]. (23) 

(Ts[no-] denotes the kinetic energy functional of the KS system.) The total ground-state 
energy Eq of the interacting system can be calculated from 

1 r ^ r ,n(r)n(r') 



^0 = T.^j,.e,,.--^ JdrJ dr'^^^^^-Y^ J dr l^xc,.(K], r)n.(r) + ExcK]. (24) 

Concerning the XC potentials we apply the L(S)DA 

^xcK] ~ J drn{r)excin+ir),n^{r)). (25) 

The most important parameterizations for the XC energy density exc{n+, nJ) (or Sxci'^^s-, 0) 
used in 2D calculations were introduced in Sec. I. 
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IV. CALCULATION OF XC POTENTIALS 



The first step of our concept to obtain XC energy densities for 2D systems is the cal- 
culation of XC potentials from GS densities. Before presenting our method we will briefly 
review the literature on the inversion of the KS equations for 3D systems. In a pioneer 
work Almbladh and Pedroza adapt parametrized XC potentials to the electron densities of 
light atoms.— An alternative approach by Aryasetiawan and Stott formulates the problem 
in terms of (A^ — 1) coupled non-linear differential equations.— Holas and March derive a 
solution by applying the Pauli potential and energy and (A^ — 1) Euler equations leading 
to a differential equation (DEQ) for the density amplitude ^Jnij)'^ Gorling describes an 
approach based on the linear response of potentials on small density modifications.-^ An 
iterative method to construct KS orbitals and XC potentials for a given electron density is 
presented by Wang and Parr^°. They use the inverted KS-Schrodinger equation to generate 
an improved effective potential for the next iteration. In another paper which mainly fo- 
cuses on kinetic energy functionals Zhao et al. propose a method which is based on Lagrange 
multipliers to gain XC potentials.— Tozer et al. train neural networka^i using XC potentials 

calculated with Zhao et al.'s method and determine fit parameters for XC functionala^i. 
□ 

Following Ref. j20|| Leeuwen and Baerends present a modified form of the inverted KS- 
Schrodinger equation's which will also be applied for a 2D system in the present paper and 
be discussed in detail below. 

The KS-Hamiltonian of axially symmetric 2D QDs reads 

+ ^rn*uy + Vn{r) + VxcAKl^) (26) 

Inversion of the KS equations means the calculation of the exact XC potentials up to a 
gauge-constant 

VxcAr) ■■= VxcAKlr) + c„ (27) 

of a system with given electron spin densities ?2o-(r), i.e. we want to find the self-consistent 
solution of the KS-Schrodinger equation under the constraint no-(r) = nKs,o-(r). The 
iteration scheme^ is based on a method presented in Ref. [2^ . It should not be considered 
as a strict proof but as plausibility argument. 

In a first step we split off the scalar XC potential from the KS-Hamiltonian ()2(i|l and 
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denote the rest by Hq^o- 

Hs,a = Ho^a- + VxC,a- (28) 

Using the eigenfunctions (y9j_o-(r) and the eigenvalues Sj o- of the KS-Hamiltonian we obtain 
the link between exact spin densities and exact XC-potentials 

VxcAr) = ^ E 7„.^,>(r) (s,,^ - i^o,.) ^,,.(r). (29) 

We denote eigenfunctions, eigenvalues, and densities of the k-th iteration step by ipj^l{r), e'j^l, 
and n^l ^.(r), i.e. they are solutions of the KS-Schrodinger equation with the KS-Hamiltonian 
-^s^o- = ^o,a + ^xcct"*- Thus, we are able to construct the XC potential for the next step 
by applying Eq. (j2l 



Vitir) = E 7..^5f (r) (^S - Ho,.) ^U^^- (30) 

With wavefunctions V^j>(r) being solutions of h'^^^^\ the iteration scheme for scalar 
potential^ is given by 

y^lAr) = Vi^-l\rf^^^. (31) 



During the iteration we will assume o- > 0- density in the {k — l)-th iteration 

step is locally too high (low) the potential of step k will be reduced (increased) at the same 
place. Consequently, the new density will be larger (smaller). However, the gauge constants 
Co- of the potentials V-^q^ cannot be determined by this method (see Ref. j^). We will 
discuss this problem in detail in the following Section. 

In the context of the iteration process the aspect of representability of the exact GS 
density is also tested: As result of a converging iteration we obtain an effective potential so 
that the KS density reproduces the density of the exact system. 

We tested the numerical procedure for a QD with six electrons and a confinement potential 
of 3meV. In the GS the two lowest shells are occupied, i.e. the system is unpolarized and the 
paramagnetic current density vanishes due to L = 0. In this test calculation we started from 
a GS density which was calculated using conventional DFT/LDA and TC parameterization 
for XC energies. Thus, we could exclude any problems arising from representability. The 
result of this test calculation shows perfect agreement between initial and final densities and 
effective potentials proving the validity of our program. 
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V. CALCULATION OF XC ENERGY DENSITIES 



In this Section we focus on the second step of the system, the extraction of XC energy 
densities from exact XC potentials. In order to establish a relation between these two 
quantities we apply the L(S)DA (|^ on the general Eq. for the functionals and obtain 

V^cAr) = excir) + n{r)^{r). (32) 

We solve for ^^{r) 

^(r) = ^ (-5xc(r) + V^cAr)) (33) 

and plug the result into the derivative of the XC energy density as a function of the radius 

^gxc(^) _ ^gxc dn^ /g4\ 
dr ^ dria dr 

After some algebra we arrive at a linear DEQ 

dsxcir) , d\ogn 

with the inhomogeneity 

I{r) = (VxcAr) + c.) ^{r) (36) 

which contains all the information about the dependence on the spin densities. In Eq. flH(i|l 
we take into account that the gauge constants Ca of the scalar potentials are not known from 
the previous step. They will be calculated later. The solution of the homogeneous part of 
DEQ (jnni) is given by ^xc^l'^) — Oi/n{r), a special solution can be calculated using the ansatz 
s^^Q{r) = /3(r)/n(r). The function (3{r) results from an elementary DEQ 

^(r)=n{r)I{r), (37) 

whose solution is 

r 

f3{r) = f dr'n{r')I{r') - (3{0). (38) 



Thus the general solution of DEQ (jH^j) is 



^xc(r) = + eTir) = + ^ / dr' n{r')I{r'). (39) 
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After substituting the inhomogeneity the result reads 

.xc(r) = ^ ^ ] V^cAr'mr'). (40) 

With the modulus of the XC energy being finite in L(S)DA 

|-ExcK]| ~ 

we choose a — P{0) — n„{0) Z^o-Co- = thus avoiding any divergent contributions. Conse- 
quently the analytical solution of (jH^j) satisfying the physical boundary conditions is 

T 

^xc(r) =T.c.^ + ^/ dr'E V.cAr'mr'). (42) 

The last step is the calculation of the gauge constants of the scalar potentials. One 
condition which has not been used yet is the agreement of the DFT GS energy and the exact 
GS energy. If we eliminate Vxc,a{^) in Eq- (1^4^ and write it in a modified form 

Egz = dTV*^{r) + K(r)) y.,.(r) + U dTVn{r)n{r) + i?xcK] (43) 

JO- \ / 

we can calculate the (exact!) XC energy Exc[na-] for the GS densities {ria-): the KS- 
wavefunctions are known from the selfconsistent solution of the KS equations (see IIV|1 . 
Thus, the expectation values / dry9*^(r) (g^^P^ + ^(^)) ^jcr{^) of the non-interacting sys- 
tem can be calculated. The Coulomb energy and the GS-energy are uniquely determined by 
the interacting system. On the other hand the XC energy for QDs in L(S)DA is given by 

oo r 

i^xcK] ~ / drn(r)£xc(r) = Ec<.A^o + / / E ^xc,.(r')^(r'), (44) 

0°" 

what makes J^a CaN^ accessible. In the case of unpolarized systems c := C| = Cj and for full 
polarization of the unoccupied spin direction is irrelevant (as the corresponding A^^^ = 0). 
Then the results for the XC energy density are unique. For partially polarized systems 
uniqueness of the results can be achieved by additionally demanding asymptotic agreement 
of Vxc,a{r) for finite systems^ 

lim Vxc,^{r) = lim VxcA^). (45) 

r— >oo " r—i-oo ^ ' ^ ' 

As a result of the inversion of the L(S)DA formalism we obtain the XC energy density 
as a function of the radius. After eliminating the radius by using the spin densities n„{r), 
or alternatively rs{r) and ^(r), we arrive at the standard representation of the XC energy 
density as a function of (r^,^). 
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J dr n{r)excir) 



drn(r)|£xc('^)| < oo 



(41) 



VI. RESULTS 



After introducing the methodology in the two previous Sections we summarize the nu- 
merical results for XC potentials and XC energy densities here. The investigated systems 
are classified due to their degree of polarization. For unpolarized = 0) and fully polarized 
systems = 1) the XC energy density is only a function of the density parameter Vg whereas 
for partially polarized systems exc depends on both and polarization ^. But first we will 
examine an analytically solvable system for two electrons to study the asymptotics of XC 
potentials. 

A. The two-electron system 

The system with two electrons will be used to find the asymptotic behavior of densities, 
KS wavefunctions, and XC potentials as QD helium can be calculated analytically for special 
values of the confinement potential. For a system with confinement energy u = 2 Ry the 
analytical expression for the singlet GS wavefunction (GS energy 6 Ry) reads 

^(ri,r2) = v^cMdri + r2|/2) V9rci(|ri - r2|) 

= CcMe-l'-^+^^l'/^ ■ aei(l + r)e-l'--^^l'/^ (46) 

with CcM = \/2/7r and Crei = 1/ \J2ti{\/2ti + 3).'^'' Consequently, we are also able to calculate 
the density, the Hartree potential and the XC potential analytically. As the full expressions 
of these quantities are very complex and do not give much insight we restrict ourselves to 
the formulas which yield the correct asymptotics for the limit of large radii. The density 
decays exponentially with 

\v^n{r) = 2'nCl^Clf'e~^' (47) 
and the Hartree potential for large r is that of a point charge 

lim l/nfr) = 4/r. (48) 
Thus, the asymptotics of the KS wavefunction R{r) = y^7rn(r) is 

lim R{r) = V27rCcMCreire-'^'/l (49) 
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R{r) = W7rn(r) can be used to calculate the XC potential 



VScc(r)/Ry = £o-r'-VH(r) + 



r dr 



(ri) R{r) 



(50) 



R{r) 



(r^ is the parabolic confinement in Ry.) The KS energy Eq = 4Ry was chosen so that the 
the XC potential vanishes in the limit r oo 



The leading term — 2/r is a manifestation of the strong selfinteraction in the case of two 
electrons. The selfinteraction potential is exactly compensated by the X potential Vx{r) = 
V}i{r)/2 in the case of a singlet configuration.— 

In summary we found from this analytical calculation that the asymptotic behavior of the 
XC potential is a rational function which converges to zero. XC energy densities cannot be 
extracted from this kind of system because of the dominating selfinteraction contribution. 

B. Unpolarized systems 

Unpolarized configurations without angular momentum are typically closed-shell systems 
with 2, 6, 12, etc. electrons. Unfortunately, the GS densities from ED for 12 and more parti- 
cles are not convergent so far, so that they cannot be used as reference densities at present. 
On the other hand, ED will give excellent results for two particles. However, predominating 
self-interaction effects in a two-electron system are a problem for the extraction of reliable 
XC energy densities. Therefore, we will focus on the six-electron system with two closed 
shells. 

To illustrate the inversion of the KS equations we investigate a six-electron system with 
a standard confinement potential 3.32 meV (see Fig. O)). The reference GS density from ED 
is depicted in (a). The quality of the ED results was tested using an increasing number of 
Slater determinants. For the present calculation with 326120 Slater determinants we found 
a good convergence of GS densities up to a radius of 8 ao. 

By applying the iteration scheme pij) we are now able to calculate the XC potential up 
to a constant (see black dotted line in (b)) which exactly reproduces the reference density 
from ED. For reasons of comparison there is also the TC XC-energy density plotted in (b) 
which was calculated with the reference density. In the next step we solve the DEQ to 



lim Vxcir) = -2/r + 1/r^ 



(51) 
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calculate the XC energy density as a function of the radius and determine the gauge constant 
Co-. The extracted XC energy density (black dashed line) is shown together with the TC 
XC energy density (grey solid) in (c). In (d) the same two curves are plotted, however, we 
eliminated the radius in favor of the density parameter. The redundancy in the extracted 
XC energy density (black, dashed curve) is explained as follows: Close to the local density 
maximum we get two exc values per density value as reflected in Fig. El (d). The deviation 
between the two values will be discussed later. 

First let us focus on the difference in the high-density region in the center and in the 
low-density regime at the edge (see Fig. IHl). While the results for high densities are reliable, 
the differences at the edge of the dot are a consequence of the finite basis set in ED. However, 
in the analytically solvable model for two electrons the XC potential (j5(Jj) converges to zero, 
whereas this is not true for the ED result (see Fig.Efb)). Consequently, the problem emerges 
again in the XC energy density (see Fig. El^c), (d)). Thus, we can reassure that the behavior 
of Vxc('^) at the edge is not due to the divergence of {dn{r) / dr) / n{r) or non-local effects. 
However, we want to emphasize that the difficulties of XC energy densities at the edge do 
not affect the accuracy of the gauge constant Ca as in Eq. £:xc('") is weighted with the 
exponentially decreasing density. 

Now let us turn to the redundant structure in sxc of Fig. El (d). As the XC energy densities 
were extracted under the assumption of LDA the order of magnitude of non-local effects is 
mapped in Fig. Efd). These problems have to be taken into account if we interpret the 
results of extracted XC energy densities. The best test for the rehability is the calculation 
of a number of systems with different confinement potentials. In Fig. |31 we summarized the 
extracted exc from systems with confinement potentials between 3 meV and 100 meV. The 
Sxc were taken from the center and the density maximum. Thus, they are not influenced by 
any edge effects and in addition they represent upper and lower boundaries for the extracted 
XC energy densities (see Fig. El (d)). All curves show the difference between a specific XC 
energy density and the TC XC energy density. While the other parameterizations from 
literature (AMGB, ISI) do not deviate too much from the TC reference both, the upper and 
lower bound of extracted XC energy densities exhibit systematic corrections in the regime 
of small density parameters or high densities. This trend is stable and consistent for many 
different strengths of confinement potentials and therefore it is reliable. The tendency to 
smaller XC energy densities is also consistent with an increasing difference of DFT and 
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ED-GS energies with growing external confinement. 

C. Fully polarized systems 

In the case of full polarization we consider the configurations with N = 6, S = 3, L = 
and iV = 8, S' = 4, L = (see Fig. |21for = 8). In analogy to the unpolarized systems we 
invert the KS equations and solve the DEQ (j35j) . The resulting XC energy densities from 
the extrema (center, maximum, minimum) are plotted in Fig. (a) for six electrons and 
in (b) for eight electrons. The strength of the confinement potential varies between 3 meV 
and 100 meV. The plots show again the deviation from the TC parameterization. As in 
the case of unpolarized systems we find for six and eight electrons systematic corrections 
in the high density regime towards smaller XC energy densities. In Fig. El (c) we check 
the consistency of the extracted XC energy densities by comparing not only results for one 
system with different confinement potentials but also different systems. The good agreement 
is an important proof for the quality and the reliability of our results. 

D. Partially polarized systems 

Finally we investigate the partially polarized systems A^ = 4, 5' = 1,L = and = 6, 
S = 2, L = 0. As an example we consider four electrons in a confinement potential of 3 meV 
(see Fig. ini). The spin-down and spin- up Fig. IHfa) densities from ED are exactly reproduced 
by the XC potentials plotted as dashed, black lines in Fig. IHl (b). Please note that the XC 
potentials for the two spin-directions exhibit the same asymptotic behavior as demanded. 
The grey, solid lines in (b) are the (gauged) TC-XC potentials as they can be calculated 
from the spin-densities in (a). After solving the DEQ (jH^I) we are able to plot the XC energy 
density as a function of the radius (see (d)). However, for partial polarization we have to 
take into account that exc depends on both, and ^. 

In order to plot the results (Fig. Ej) we have to collect all combinations {rs,exc) for given 
polarization ^ from the systems with four and six electrons and confinement energies from 
3 meV to 100 meV. All the values were taken from areas with a reasonably high density. Only 
in this regime results are reliable as follows from the same arguments as for unpolarized and 
fully polarized systems. The summary of extracted XC energy densities is shown in Fig. [7| 
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for the polarizations ^ = 0.2 (a), ^ = 0.4 (b), ^ = 0.6 (c), and ^ = 0.8 (d) as deviation 
from the TC parameterization. Similar to unpolarized and fully polarized systems we find 
systematic corrections for small compared to TC, AMGB, and ISI. The reliability of the 
extracted exc is emphasized by the fact that the difference between them is not larger than 
for ^ = and ^ = 1. 

VII. SUMMARY AND OUTLOOK 

In this paper we present a new approach how to extract 2D-XC energy densities from 
GS densities and energies of QDs calculated by ED. Our focus was on configurations with 
arbitrary polarization but without paramagnetic current density. In comparison with pa- 
rameterizations from literature (TC^"^, AMGB^, ISI^I) we find systematic corrections in the 
regime of small density parameters or high densities for all polarizations. This result was 
consistently accomplished by exploiting the redundancy arising from different configurations 
and/or different strengths of the external confinement potentials. 

However, our results are restricted to a small range of density parameters. Therefore, it 
would be useful to extend the calculations to more confinement potentials or to confirm them 
by an alternative method which also provides an approach for an improved parameterization. 

Another open question is the role of the corrections due to the paramagnetic current 
density or vorticity^>^ which lead to an XC vector potential. Up to now these effects have 
been considered in form of complicated interpolations.— >2^^ This problem will be tackled 
in a forthcoming paper. 

VIII. ACKNOWLEDGEMENT 

A. W. thanks the the RRZE Erlangen and the DFG (Ro 522/19-1). 



1 P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 (1964). 

2 W. Kohn and L. J. Sham, Phys. Rev. 140, A 1133 (1965). 

3 U. von Barth and L. Hedin, J. Phys. C 5, 1629 (1972). 

^ U. K. Rajagopal and J. Callaway, Phys. Rev. B 7, 1912 (1973). 



15 



^ G. Ortiz, D. M. Ceperley, and R. M. Martin, Phys. Rev. Lett. 71, 2777 (1993). 

6 M. Koskinen, M. Manninen, and S. M. Reimann, Phys. Rev. Lett. 79, 1389-1392 (1997). 

^ D. G. Austing, S. Sasaki, S. Tarucha, S. M. Reimann, M. Koskinen, and M. Manninen, Phys. 

Rev. B 60, 11514 (1999). 

^ K. Hirose and N.S. Wingreen, Phys. Rev. B 59, 4604 (1999). 

9 O. StefFens and M. Suhrke, Phys. Rev. Lett. 82, 3891 (1999). 

1° A. Wensauer, O. StefFens, M. Suhrke, and U. Rossler, Phys. Rev. B 62, 2605 (2000). 

1^ A. Wensauer, J. Kainz, M. Suhrke, and U. Rossler, phys. stat. sol. (b) 224, 675 (2001). 

M. Ciorga, A. Wensauer, M. Pioro-Ladriere, M. Korkusinski, J. Kyriakidis, A. S. Sachrajda, 

and P. Hawrylak, Phys. Rev. Lett. 88, 256804 (2002). 

13 A. Wensauer, M. Korkusinski, and P. Hawrylak, Phys. Rev. B 67, 035325 (2003). 

1^ B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 (1989). 

1^ C. Attaccalite, S. Moroni, P. Gori-Giorgi, and G. B. Bachelet, Phys. Rev. Lett. 88, 256601 
(2002). 

16 C.-O. Almbladh and A. C. Pedroza, Phys. Rev. A 29, 2322, (1984). 

1^ F. Aryasetiawan and M. J. Stott, Phys. Rev. B 38, 2974 (1988). 

18 A. Holas and N. H. March, Phys. Rev. A 44, 5521 (1991). 

19 A. GSrIing, Phys. Rev. A 46, 3753 (1992). 

20 Y. Wang and R. G. Parr, Phys. Rev. A 47, 1591 (1993). 

21 Q. Zhao and R. C. Morrison and R. G. Parr, Phys. Rev. A 50, 2138 (1994). 

22 R. van Leeuwen and E. J. Baerends, Phys. Rev. A 49, 2421 (1994). 

23 D. J. Tozer, V. E. Ingamells, and N. C. Handy, J. Chem. Phys. 105, 9200 (1996). 

24 D. J. Tozcr, N. C. Handy, and W. H. Green, Chem. Phys. Lett. 273, 183 (1997). 

25 A. K. Rajagopal and J. C. Kimbah, Phys. Rev. B 15, 2819 (1977). 

26 J. p. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
2''' M. Seidl, University of Regensburg, unpublished (2001). 

28 S. M. Reimann, M. Koskinen, and M. Manninen, Phys. Rev. B 62, 8108 (2000). 

29 A. Wensauer, M. Korkusinski, and P. Hawrylak, submitted to Phys. Rev. B, 2003. 

30 O. Steffens, U. Rossler, and M. Suhrke, Europhys. Lett. 42, 529 (1998). 

31 K. Capelle and G. Vignale, Phys. Rev. Lett. 86, 5546 (2001). 

32 G. Vignale and M. Rasolt, Phys. Rev. Lett. 59, 2360 (1987). 

16 



G. Vignale and M. Rasolt, Phys. Rev. B 37, 10685 (1988). 

D. Levesque, J. J. Weis, and A. H. MacDonald, Phys. Rev. B 30, 1056 (1984). 

G. Fano and F. Ortolani, Phys. Rev. B 37, 8179 (1988). 

M. Rasolt and F. Pcrrot, Phys. Rev. Lett. 69, 2563 (1992). 

M. Taut, J. Phys. A 27, 1045 (1994). 

C. Fihppi, C. J. Umrigar, and M. Taut, J. Chem. Phys. 100, 1290 (1994). 



17 



FIG. 1: Comparison of TC™, AMGB^, and ISI^l parameterizations as a function of density 
parameter and polarization. In (a) and (b) the relative deviation 

4c ^^/^xc - 1 between TC and 

AMGB is plotted on different scales for r^. (c) and (d) show the difference between ISI and TC 
(eAMGB/^TC _ 1) and AMGB (^^^^/e^g - 1). 

FIG. 2: Comparison of GS densities from ED and DFT for a fully polarized eight-electron system 
(L = 0, S" = 4) with pronounced differences in the density extrema. The three DFT curves with 
TC, AMGB, and ISI parameterization lie so close that they cannot be resolved. 



FIG. 3: Numerical results of an unpolarized six-electron system with L = S = and confinement 
potential 3.32 meV. In (a) we show the GS density from ED in (a) as a function of the radius. 
It is reproduced by the XC potential (black, dashed line) in (b) which gives rise to the extracted 
XC energy density (black, dashed line) in (c). In (d) the extracted XC energy density is plotted 
versus the density parameter (black, dashed line). The redundant structure is a conseuquence 
of the density profile (see text). The grey, solid line in (b), (c), and (d) is the corresponding result 
calculated with TC parameterization. 
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FIG. 4: Summary of extracted XC energy densities for unpolarized systems (six electrons, different 
confinement potentials) as a function of r^. We plotted the deviation of the AMGB, the ISI, and 
the extracted XC energy densities from the TC parameterization. The results from the density 
maximum and the dot center are upper and lower bounds for the calculated exc-values. For small 
Ts we find a systematic correction towards smaller XC energy densities. 

FIG. 5: Summary of extracted XC energy densities for polarized systems (six (a) and eight (b) 
electrons, different confinement potentials) as a function of Vg. We plotted the deviation of the 
AMGB, the ISI, and the extracted XC energy densities from the TC parameterization. The results 
from the density extrema represent upper and lower bounds for the calculated £xc-values. (c) 
shows the combined results for six and eight electrons and different confinement strength. For 
both electron numbers, we find consistent and systematic corrections towards smaller XC energy 
densities at small density parameters. 

FIG. 7: Summary of extracted XC energy densities for partially polarized systems (combined results 
for four and six electrons, different confinement potentials) as a function of and polarization 
^ = 0.2 (a), 0.4 (b), 0.6 (c) and 0.8 (d). We plotted the deviation of the AMGB, the ISI, and the 

extracted XC energy densities from the TC parameterization. For all systems we find consistent 
and systematic corrections towards smaller XC energy densities at small density parameters. 



FIG. 6: Numerical results of a partially polarized four-electron system with L = 0, 5 = 1 and 
confinement potential 3 meV. In (a) and (b) we show the GS densities from ED for up- and down- 
spins as a function of the radius. They are reproduced by the XC potentials (black, dashed lines) 
in (c) which give rise to the extracted XC energy density (black, dashed line) in (d). The grey, 
solid line in (b), (c), and (d) is the corresponding result calculated with TC parameterization. 
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